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ITA) JPEG 压缩 域 脆弱 图 像 水 印 算法 
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摘 要 : 为 了 解决 JPEG 图 像 水 印 算 法 在 上 座 入 容量 、 不 可 感知 性 、 实 时 性 和 窜改 定位 的 问题 ， 提 出 了 一 种 基于 炳 编码 
的 JPEG 压缩 域 脆弱 图 像 水 印 算 法 。 首 先 选 择 JPEG J8 2822 MRT KERA, RAER UD SEE fe GE 9) dfeiÉ a 89 DCT 和 量 
化 运算 ， 保 证 算法 在 谨 入 阶段 对 图 像 的 改变 量 较 小 且 有 具有 较 好 的 实时 性 ; 接着 根据 业 编 码 过 程 中 某 些 比特 不 会 参与 哈 
夫 曼 编码 这 一 原则 ， 在 其 中 选择 符合 条 件 的 系数 诅 入 水 印 ， 将 水 印信 息 砍 入 到 这 些 系 数 的 最 低 有 效 位 中 ， 以 规避 修改 
哈 夫 曼 编 码 系 数 后 带 来 的 编码 误差 ,进一步 加 强 了 算法 的 黎明 性 。 实 验 表 明 ， 该 算法 不 仅 具 有 较 大 的 详 入 容量 和 很 好 
的 不 可 感知 性 ， 而 且 脆 弱 性 较 高 且 能 够 对 恶意 窜改 进行 准确 定位 。 
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JPEG compressed domain fragile image watermarking algorithm based on entropy coding 
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uu Abstract: In order to solve the problem of embedding capacity, imperceptibility, real-time and tampering localization in JPEG 


image watermarking algorithm, this paper proposes a JPEG watermarking algorithm based on entropy coding. Firstly, the 


algorithm locates the process of JPEG entropy coding to embed the watermarking, which can avoid the forward and reverse 


quantization operations and ensure that during embedding image changes is small and has good real-time performance. Secondly, 


according to the principle that some bits in the entropy coding process do not participate in the Huffman coding, the algorithm 
chooses appropriate coefficients to embed watermark. Our algorithm embeds watermarking information into the least significant 
bits of these coefficients to avoid the coding error caused by the modification of the Huffman coding coefficient, which should 
further enhance the imperceptibility of the algorithm. Experiments results show that the algorithm not only has a large embedded 
capacity and good imperceptibility, but also has high fragility and can accurately locate malicious tampering. 
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0 引言 JPEG 脆弱 水 印 算法 ,根据 8X8 图 像 块 中 的 DCT 系数 生成 4 个 
比特 的 水 印信 息 ， 然 后 根据 密 钥 随 机 的 众 入 到 相 邻 的 4 个 图 像 

随 着 互联 网 和 多 媒体 技术 的 快速 发 展 ，JPEG 格式 的 区 块 DCT 系数 的 最 低位 。 上 述 文献 中 提出 的 脆弱 水 印 算法 都 是 
得 到 了 广泛 的 应 用 ， 但 是 目前 对 于 图 像 水 印 算法 的 研究 绝 大 多 Æ DCT 系数 上 嵌入 水 印 ， 每 次 水 印 的 嵌入 和 提取 都 会 涉及 到 
数 都 是 针对 于 BMP 格式 图 像 , 针对 JPEG 压缩 域 的 图 像 水 印 算 E 向 和 逆向 的 DCT 和 量化 运算 , 算法 的 时 间 复 杂 度 较 高 。 另 外 
法 研究 较 少 。 而 其 中 大 部 分 的 研究 主要 集中 于 JPEG 压缩 域 的 虽然 上 述 算 法 对 于 常规 攻击 具有 一 定 的 脆弱 性 ， 但 是 对 图 像 质 
和 鲁 棒 性 水 印 算法 ,用 于 图 像 内 容 的 完整 性 以 及 真实 性 验证 — 量 影响 较 大 ， 不 可 感知 性 较 差 。 文 献 [12] 提 出 了 一 种 灵活 的 
NHJ JPEG 压缩 域 脆弱 性 水 印 算法 的 研究 相对 更 少 , 与 鲁 棒 性 水 JPEG 压缩 域 水 印 算法 , 通过 更 改 Zig-zag 排序 后 一 维 数组 中 最 
印 算法 相 比 ， 脆 弱 水 印 算法 在 图 像 受到 攻击 时 会 比较 敏感 。 ”后 一 个 符合 条 件 的 系数 的 位 置 来 嵌入 水 印 ， 该 算法 有 一 定 的 脆 
文献 [10] 提 出 一 种 可 逆 的 JPEG 脆弱 水 印 算法 , 通过 折 针 量化 后 。 弱 性 同时 能 抵抗 特定 的 噪声 攻击 ， 但 是 该 算法 不 能 检测 具体 的 
系数 的 哈 希 值 来 生成 水 印 ， 然 后 将 2 个 比特 分 别 骨 入 到 每 个 块 窜改 区 域 。 同 时 由 于 Zig-zag 中 的 顺序 对 于 炉 编 码 过 程 非常 重 
中 2 个 系数 的 最 低 有 效 位 中 。 文 献 [11] 提 出 基于 领域 比较 的 。” 要 ,调整 系数 顺序 间接 改变 最 后 图 像 的 压缩 结果 ， 对 图 像 质 量 
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的 影响 较 大 。 文献 [13] 提 出 了 一 种 零 失 真 的 JPEG 水 印 算法 , 通 
过 同时 修改 哈 夫 曼 码 表 和 比特 流 来 嵌入 水 印 ， 由 于 并 没有 修改 
任何 JPEG 量化 之 后 的 系数 ， 所 以 该 算法 不 会 对 图 像 造成 任何 
质量 影响 , 虽然 该 算法 脆弱 性 较 好 , 但 是 也 不 能 进行 窜改 定位 。 
综 上 ， 由 于 焙 编 码 阶段 在 量化 过 程 之 后 ， 避 免 了 正 向 和 逆 
向 的 DCT 和 量化 运算 ,能 具有 很 好 的 透明 性 和 实时 性 , 本 文 提 
Hr PhFD JPEG 脆弱 图 像 水 印 算法 。 算 法 思想 是 
首先 定位 JPEG Jis fr Be REA Dur A 5 S RE SIE EG 
特 ， 然 后 从 中 选择 交流 系数 将 水 印 典 入 到 其 最 低 有 效 位 中 。 
于 交流 系数 具有 比较 敏感 的 特点 ， 所 以 本 文 算法 不 仅 对 常见 攻 
击 具 有 良好 的 脆弱 性 ， 而 且 能 够 对 恶意 窜改 进行 准确 定位 。 


1 JPEG 编码 过 程 及 骨 入 位 置 选择 


JPEG 图 像 压 缩 算法 具体 流程 如 图 1 所 示 ， 对 于 BMP 图 像 
首先 进行 颜色 空间 转换 ， 然 后 对 每 个 颜色 分 量 进 行 8x8 分 块 ， 
接着 对 每 个 块 进行 DCT 和 量化 运算 ， 最 后 经 过 Zig-zag 排序 和 
Jj i 3 Je: LA JPEG 格式 保存 得 到 最 终 的 JPEG 图 像 。 


RGB 颜色 空间 YUV 颜 色 空 间 8*8 分 块 
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JPEG 图 像 编码 过 程 
JPEG [ERI A RA RU A A IET EE, RARO 
包括 行程 编码 、 比 特 编码 和 哈 夫 曼 编 码 这 三 个 过 程 史 。 为 了 验 


| Chinax iv (ERBTI. 
300 : vu E eH 法 


李 ROT 


的 是 每 一 个 处 理 单元 最 多 只 能 有 16 个 0, 也 就 是 说 如 果 有 超过 
16 个 0 的 情况 出 现 ， 就 需要 单独 进行 分 组 。 如 果 最 后 一 个 单元 
都 是 0， 那么 就 使 用 一 个 特殊 的 标识 符 EOB 来 表示 。 行 程 编码 
过 程 如 图 3 所 示 。 
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ORANE 35,7,0,0,0,-6,-2,0,0,-9,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,8,0,0,0,...,0 
35 7 0,0,0,-6 00.9 | 0,0,..,0,8 0,0,...,0 
ememm | 35 | 7 | ooo | | 00.9 | 00，0 008 | 00 0 
039 | (07) | G9 | (orz2) | e» | (5,0) 28) | EoB 


图 3 行程 编码 过 程 

行程 编码 又 具体 分 为 3 步 ， 第 一 步 是 将 非 0 和 该 数 前 面 的 
0 划分 为 一 个 处 理 单元 ,第 二 步 是 对 每 一 个 单元 中 有 超过 15 个 
0 的 情况 进一步 分 组 ， 第 三 步 对 以 上 单元 再 进行 编码 得 到 一 个 
二 元 组 。 二 元 组 左边 的 参数 代表 0 的 个 数 ， 右 边 的 参数 代表 非 
0 系数 的 值 。 例 如 (3,-6) 代表 -6 前 面 还 有 3 个 0，(15,0) 表示 
0 前 面 还 有 15 个 0， 所 以 (15,0) 表示 该 分 组 有 16 个 0. 
在 行程 编码 结束 后 ， 将 首先 对 每 个 二 元 组 的 右边 参数 进行 
比特 编码 , JPEG 对 此 提供 了 一 张 标准 的 码 表 05。 这 里 比特 编码 
的 本 质 就 是 将 目前 的 十 进 制 数 转换 为 其 二 进 制 编码 并 记录 其 二 
进 制 码 长 度 ， 如 果 参 数 为 正 数 ， 那 么 编码 结果 就 是 其 二 进 制 编 
码 原 码 ; 如 果 参 数 是 负数 , 那么 就 是 其 正 数 二 进 制 编码 的 反 码 。 
比如 行程 编码 后 得 到 二 元 组 (3,-6), 这 时 先 对 -6 进行 比特 编码 ， 
-6 为 负数 所 以 其 编码 结果 为 其 正 数 的 二 进 制 编码 的 反 码 ， 其 正 
数 6 的 三 进 制 编码 为 110， 所 以 -6 的 比特 编码 为 001， 长 度 是 
3， 那 么 比特 编码 第 一 步 的 结果 就 为 (3,3,001)。 其 中 三 元 组 
(3,3,001) 中 参数 从 左 到 右 分 别 代 表 0 的 个 数 、 行 程 编码 后 二 
元 组 右边 参数 二 进 制 编码 的 长 度 以 及 甚 二进制 编码 。 这 里 需要 
注意 的 是 ,在 行程 编码 中 专门 说 明 每 一 个 分 组 内 最 多 仅 有 16 个 
0， 同 时 用 〈15,0) 代表 连续 出 现 16 次 0 的 情况 ， 所 以 三 元 组 
最 左边 的 参数 不 会 超过 15。 并 且 JPEG 编码 规则 中 规定 ， 比 特 
编码 第 一 步 处 理 的 右边 参数 只 会 介 于 -2047=-21 到 204722!! 之 


证 散 入 位 置 的 显著 性 及 透明 性 ， 必 须 先 介绍 JPEG 图 像 压缩 血 
EBD. EALAR DCT, EAKR Zig-zag 
HEFE, KPIREN EEES AKE 64 的 一 维 
数组 ， 如 图 2 所 示 。 


ES, 所 以 其 二 进 制 位 数 最 多 为 11 位 , 也 就 意味 着 三 元 组 中 间 
的 参数 不 会 超过 11。 那么 这 时 可 以 将 三 元 组 的 前 两 个 部 分 合 
起 来 ， 用 一 个 字 节 来 表示 ， 即 高 4 位 表示 出 现 0 的 次 数 ， 低 4 
位 表示 编码 位 数 。 比 特 编码 过 程 如 图 4 所 示 。 


加 原始 数据 35,7,0,0,0,-6,-2,0,0,-9,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,8,0,0,0,...,0 


图 2 Xi 8s cd 

可 以 看 到 每 个 8x8 的 分 块 在 经 过 DCT 运算 尤其 是 量化 之 
后 ,和 矩阵 中 出 现 0 的 次 数 是 非常 多 的 。 由 于 一 个 分 块 在 经 过 DCT 
变化 后 ， 其 直流 系数 会 代表 其 原 有 分 块 的 大 部 分 能 量 ， 一 般 值 
都 比较 大 ， 而 右 下 角 剩 余 的 交流 系数 则 代表 了 原 有 分 块 的 细节 
分 量 ， 一 般 值 都 比较 小 ， 且 越 靠 右 下 角 值 越 小 ， 所 以 再 经 过 量 
化 后 矩阵 右 下 角 区 域 基本 都 为 0。 那 么 接着 将 这 个 分 块 中 量化 
后 的 系数 按照 Zig-zag 排序 后 形成 的 一 维 数组 中 , 0 的 位 置 会 相 
对 比较 集中 ， 故 焙 编 码 过 程 的 第 一 步 就 是 对 这 些 0 进行 处 理 。 
为 了 节省 编码 位 数 ，JPEG 使 用 了 行程 编码 (RLE) 将 数组 中 的 
非 0 和 该 数 前 面 的 0 的 个 数 作 为 一 个 处 理 单元 。 这 里 需要 说 明 


m 


| 四 原始 政 据 | 35,7,0,0,0,-6,-2,0,0,-9,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,8,0,0,0,...,0 
| 35 [| 7 [oos] -2 oo. | 0,0,0,8 [oo_o 
四 RLE 编 码 | 35 | | 0,0,0,-6 | 2 0,0,-9 | 0,0,...,0 | 0,0,8 | 0,0,...,0 
| (0,35) ^| (07) | (3,-6) T (0-2) i 2,-9) mI (15,0) | (2,8) | EOB 
(0,6, 100011) | (0,3, 111) | (3,3, 001) Toz , 01) (2,4, 0110 ! 15,-) ] 2,4, 1000) T EOB 


(0x 6,100011) | (0x3,111) | 0x33,001) | (0x2,01) | (0x24,0110) | (OxF0,-) | (0x24,1000) | EOB 


图 4 比特 编码 过 程 

接 下 来 还 需要 进行 哈 夫 曼 编 码 ， 即 还 需要 对 比特 编码 后 得 
Sil — 76H B 265 21 CIE TT I DER AR, iX HE REI A A 
过 程 的 最 后 一 步 。 JPEG 提供 了 两 张 标准 的 码 表 , 分 别针 对 于 直 
流 系数 和 交流 系数 , 具体 码 表 可 参见 文献 [16], SEA AI ETE 
如 图 5 所 示 。 
丸 为 直流 分 量 比 交流 分 量 更 加 重要 ， 将 水 印 嵌 入 到 量化 后 
的 直流 分 量 当 中 会 对 载体 图 像 的 质量 产生 较 大 影响 0% 4， 所 以 
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本 文选 择 将 水 印信 息 柑 入 交流 系数 当中 。 在 图 5 中 即 从 系数 7 
开始 即 可 算 作 交 流 系 数 , 在 比特 编码 时 间 其 编码 值 为 (0x3,111)， 
从 编码 过 程 图 5 可 以 看 出 ， 比 特 编 码 右边 系数 编码 值 111 不 会 
参与 哈 夫 曼 编 码 ， 所 以 修改 右边 参数 的 值 不 会 造成 哈 夫 曼 码 纺 
码 的 错位 ， 这 样 嵌 入 水 印 后 的 图 像 将 具有 很 好 的 透明 性 ;同时 
这 种 修改 也 不 涉及 正 向 或 者 逆向 的 DCT 和 量化 操作 ， 所 以 算 
法 时 间 具 有 很 好 的 实时 性 。 并 且 由 于 每 个 分 块 都 会 嵌入 水 印 ， 
因此 在 受到 窜改 攻击 时 仅 会 影响 被 窜改 的 块 ， 而 未 受 窜改 的 块 
不 会 受到 影响 ， 所 以 可 以 依据 分 块 准确 定位 窜改 。 综 上 ， 本 算 
法 选择 每 个 块 中 交流 系数 在 比特 编码 时 期 二 元 组 的 最 右边 参数 
作为 水 印 的 嵌入 位 置 。 


(eet 35,7,0,0,0,-6,-2,0,0,-9,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,8,0,0,0,...,0 
35 7 0,0,0,-6 -2 0,0,-9 0,0,..,0,8 0,0, 
,0 
J@RLE 编 码 | 35 7 0,0,0,-6 2 0,0,-9 0,0, 0,0,8 0,0, 
Lj 0 
0,35) (0,7 3,-6) 0,-2) (2,-9) a50) Q8 EOB 
(0,5, 100011) | (0,3, 111)| (3,3, 001) (0,2 (24,0110) | (15,- 4, 1000) | E08 
| Gema 01) 
(0x6,100011) | (0x3,111 (0x33,001) |(0x2,01)| (0x24,0110) | oxFo | (0x24,1000) | E08 
- 
| 100 | 100011 | 100 |111| 1111 |001|01 | 01 | 1111 |0110| 1111 | 1111 | 1000 | 1010 
angge 
B 1111 1111 |1111 
0101 0100 001 |0100 
c 
1001 000111001111111 111110101001010111111111010001101111111 100111111111010010001010 
全 序列 化 ——— E —— 1 


91 CF FE AS 7F D1 BF CF FA 45 | 
kp 5 Ji ERE 

于 本 算法 是 在 编码 过 程 嵌 入 水 印 ， 需 要 在 解码 过 程 提 取 
水 印 ， 那 么 即使 是 脆弱 水 印 算法 也 必须 对 编 解 码 过 程 具 有 和 鲁 棒 
性 ， 如 果 选 择 最 后 一 个 交流 系数 嵌入 水 印 ， 该 系数 过 于 脆弱 ， 
算法 在 解码 部 分 就 无 法 正确 提取 水 印 , 不 能 达到 水 印 提取 要 求 。 
故 本 算法 选择 第 一 个 交流 系数 嵌入 水 印 ， 并 且 在 每 个 块 中 嵌入 
次 来 达到 对 编 解码 过 程 的 鲁 棒 。 


2 ”水 印信 息 的 生成 


为 了 进行 水 印 嵌 入 ， 首 先 需要 根据 水 印 图 片 生 成 可 用 的 水 
印 序列 。 处 理 的 操作 过 程 为 水 印信 息 的 二 值 化 、 水 印 置 乱 处 理 
和 水 印信 息 降 维 , 最终 的 目的 是 得 到 一 个 一 维 的 二 值 水 印 数组 。 
具体 步 又 如 下 : 

a) 二 值 化 ， 就 是 将 水 印 图 像 的 像素 点 值 转化 为 0 或 者 255 
的 过 程 ， 实 际 编码 过 程 中 会 映射 到 0 和 1。 首先 水 印 图 像 为 灰 
度 图 像 ， 如 果 是 彩色 图 像 需要 先 转换 为 灰 度 图 像 。 本 文中 映射 
过 程 通过 阔 值 来 进行 ， 选 取 作 为 128 为 阔 值 ， 大 于 等 于 128 的 
像素 点 值 都 会 被 映射 为 1， 小 于 128 的 像素 点 值 都 映射 为 0。 

b) Armold 变换 。Arnold 变换 是 水 印 算法 中 使 用 最 广 的 
种 置 乱 算 法 ， 本 质 是 将 数字 图 像 中 的 像素 点 重新 排列 ， 以 此 来 
打 乱 图 像 的 显示 。 本 文 使 用 的 Arnold 置 乱 具体 计算 公式 如 下 : 


x| d lx dN D 
317 bra] Pee PS ( 


其 中 : (x，y) 表示 原始 图 像 的 像素 点 ，(x'，y') 为 变换 后 的 像素 
点 ，N 为 图 像 的 高 度 

c) 水 印 降 维 , 即将 二 维 的 水 印 图 像 信 息 和 矩阵 转换 为 一 维 的 
数组 。 可 以 通过 行 转换 的 方式 ， 也 可 以 通过 列 转 换 的 方式 ， 需 
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要 注意 的 是 按 哪 些 方式 降 维 ， 今 后 需要 同样 的 方式 升 维 。 首 先 


需要 计算 出 一 维 数组 的 长 度 ， 计 算 公式 如 下 ， 


C-MxN 


其 中 : M 代表 水 印 图 像 的 高 度 , N 代表 水 印 图 像 的 宽度 。 最终， 


水 印信 息 可 以 表示 为 


W -(WG), 1<i<M xN)} 


3 KEPARA MERE A 


3.1 zKED BRA iE fe 


(3) 


本 文 基 于 透明 性 和 实时 性 要 求 选 择 在 DCT 系数 量化 及 
Zig-zag 排序 之 后 的 比特 流 焙 编码 的 过 程 中 嵌入 水 印 。 另 外 如 果 


载体 图 像 是 彩色 图 像 ， 如 果 同 时 对 Y, Cr, Cb 通 


iB BEN ZKED ADT 


图 像 质量 影响 较 大 。 由 于 JPEG 一 个 编码 单元 中 各 颜色 分 量 块 


数 的 比例 为 ，Y:Cr:Cb=4:1:1， 选 Cr Cb 通道 块 过 少 会 很 可 能 不 


WX. 


PERAREHER, MON TOR GLELR HUGE Y 分 量 进行 水 印 


了 避免 水 印 在 传输 过 程 中 可 能 发 生 错误 的 问题 ， 本 文 需 


要 对 能 嵌入 水 印 的 块 进行 筛选 .从 每 个 块 第 一 个 


交流 系数 开始 ， 


如 果 系 数 非 0 且 绝 对 值 大 于 1， 那 么 进行 标记 然后 继续 往 下 饥 


历 ， 直 到 找到 大 个 满足 条 件 的 系数 后 ， 本 文才 选择 该 块 插入 水 


印 。 柑 入 时 将 在 每 一 个 块 中 嵌入 一 个 比特 的 水 E 


水 印 比特 信息 将 被 重复 嵌入 该 块 的 磊 个 符合 条 件 的 交流 系数 中 。 


信息， 每 一 个 


水 印 的 肉 入 过 程 在 编码 阶段 完成 ， 水 印 嵌 入 的 


i 程 图 如 图 6 所 


示 。 

输入 原始 

bmp 图 像 输入 水 印 图 像 
* * 

颜色 空间 转换 并 一 

8*8 分 块 值 化 

+ 
DCT 并 量化 Arnold 变 换 
* * 
Zig-Zag 排 序 并 行程 编码 降 维 
* 
— 获得 比特 编码 后 的 系数 C,~C,， 


^| ”从 C, 开 始 ， 寻找 满 足 C+0 且 |C|>1 的 系数 C, 


能 入 水 印 到 该 系数 的 最 低 有 效 位 中 


Y 
N 
Y 


I N — ms > 


选择 下 一 个 8*8 块 


p ---X- 
哈 夫 曼 编 码 
i 
1838 d X k ENE HO JPEG [EL t 
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人体 步骤 如 下 : 
a) 读 取 原始 的 JPEG 图 像 ， 对 于 每 一 个 8x8 BER, FER 
码 阶段 ， 执 行 哈 夫 曼 解 码 和 行程 解码 以 获得 比特 编码 时 期 的 最 
ARAE, HEW ZG), ， 其 中 0<j 和 63, i=1,2...， N 。 
b) 对 于 每 一 个 8x8 的 块 , 则 需要 找到 从 第 一 个 交流 系数 开 
台 ， 满 足 非 0 且 绝 对 值 大 于 1 的 系数 ， 定 位 其 在 二 进 制 码 流 中 
的 位 置 ， 具 体 定义 如 下 所 示 : 


Pos,-Min(j) :Z(j)) £0 |Z( [>1, 1€j€63,i212,....N (4) 


此 时 ， 第 一 


个 满足 条 件 的 系数 (First Appropriate Data, FAD ) 
即 为 
FAD = Z (Pos,) (5) 

c) 从 FAD 的 位 置 开 始 依次 往 下 遍历 ， 如 果 还 有 系数 满足 
非 0 且 绝对 值 大 于 1 的 条 件 ,那么 进行 标记 然后 继续 往 下 遍历 ， 
直到 找到 个 满足 条 件 的 系数 后 ， 跳 转 到 下 如 果 在 
该 块 中 没有 找到 磊 个 满足 条 件 的 系数 ， 则 跳 过 这 个 块 ， 对 下 一 
个 块 从 步骤 2 开始 继续 进行 遍历 选择 。 

d) 对 于 每 个 选 定 的 8x8 块 ， 根 据 水 印信 息 W 的 值 来 依次 


下 一 个 步骤; 


修改 大 个 系数 的 值 ， 具 体 每 个 系数 的 修改 过 程 如 下 所 示 ; 
Ci» CA 为 偶数 并 且 W, =0 

c= C. +1, C 为 偶数 并 且 W, =1 T. 
"ens C^ 为 奇数 并 且 W -1 
C5 -1， — CK 为 奇数 并 且 W, =0 


其 中 : CiS 代表 当前 需要 嵌入 水 印 的 原 系数 值 ，u, v 下 标 代表 其 
在 矩阵 块 中 的 位 置 , AC 代表 其 为 交流 系数 。C%' 为 对 应 第 i 个 
系数 嵌入 水 印 比特 后 的 新 系数 值 ， 折 为 第 i 个 块 对 应 的 一 维 水 
印 的 比特 值 。 

e) 对 于 每 个 符合 条 件 的 8x8 Ee, 在 码 流 中 用 得 到 的 新 的 值 
CiS 依次 替换 原始 的 Co. ， 然 后 再 进行 哈 夫 曼 编 码 ， 从 而 得 到 
了 最 终 的 供 入 水 印信 息 的 JPEG 图 像 序列 。 
3.2 水 印 提取 过 程 
在 水 印 提取 阶段 ， 需 要 在 JPEG 的 解码 过 程 中 提取 水 印 。 
同样 ， 对 于 每 个 8x8 的 图 像 块 ， 都 需要 进行 哈 夫 曼 解 码 ， 然 后 
首先 定位 到 每 一 个 符合 条 件 的 块 中 的 第 一 个 非 0 且 绝 对 值 大 于 
1 的 交流 系数 (First Appropriate Data, FAD), 根据 FAD 的 位 置 
依次 往 后 遍历 ， 按 序 找到 符合 条 件 的 大 个 系数 ， 然 后 从 每 一 个 
系数 中 提取 一 个 比特 的 水 印信 息 ， 统 计 从 大 个 系数 中 提取 出 的 
水 印信 息 的 值 ， 根 据 1 出 现 的 次 数 来 决定 当前 块 中 提取 的 水 印 
信息 的 值 。 将 所 有 块 中 的 水 印信 息 提 取 后 ， 再 进行 一 维 水 印 数 
组 的 升 维 、 反 置 乱 处 理 、 以 图 片 格式 保存 ， 即 得 到 了 提取 出 的 
水 印 图 像 。 水 印 提取 的 流程 图 如 图 7 所 示 ， 有 具体 步骤 如 下 : 
a) ERRAK K JPEG 图 像 ， 对 于 每 一 个 8x8 的 块 ， 
在 粒 编 码 阶 段 ， 执 行 哈 夫 曼 解 码 到 比特 编码 以 获得 比特 编码 时 


IT 
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期 的 最 右 系 数 的 值 ， 其 值 记 为 : 
0< j<63, i=1,2,..., N., 

bo ERAZ O 中 的 第 一 个 非 0 且 绝 对 值 大 于 1 的 交流 系 
数 ， 即 FAD'. 

c) 从 FAD' 的 位 置 开始 依次 往 下 遍历 ,如 果 系 数 CL ' 满足 
非 0 且 绝对 值 大 于 1 的 条 件 , 那么 进行 标记 然后 继续 往 下 遍历 ， 
直到 找到 磊 个 满足 条 件 的 系数 后 ， 则 选中 该 块 进行 水 印 提取 ， 
进行 下 一 个 步骤 ; 如 果 在 该 块 中 没有 找到 让 个 满足 条 件 的 系数 ， 
就 跳 过 这 个 块 ， 对 下 一 个 块 从 步骤 2 开始 继续 进行 遍历 选择 。 

d) 在 被 选中 的 块 中 ， 根 据 系数 C ' 的 值 来 从 每 个 系数 
取 一 个 比特 的 水 印信 息 ， 统 计 天 个 系数 中 提出 的 工 的 次 数 ， 具 
体 公 式 如 下 : 


ZU ,. Xm. 


Gu 


Count CC ' 为 偶数 
Count= 7 n CI) 
Count+ 1 C^ 为 奇数 


其 中 : Count 初始 值 为 0， 如 果 当 前 系数 值 为 偶数 ，Count 
EDE: 如 果 当 前 系数 值 为 奇数 ，Count 值 自 增 1。 即 统计 大 个 
系数 中 1 出 现 的 次 数 ， 然 后 根据 当前 块 中 1 出 现 的 次 数 获得 当 


前 块 中 提取 的 水 印信 息 。 有 具体 公式 如 下 : 
0, Count < £ 
W'= E (8) 
l, Count 2 一 
2 


ies Count 代表 该 中 从 个 系数 中 提取 的 信息 中 1 t 


现 的 次 数 。W' 表示 从 符合 条 件 的 块 中 提出 的 第 i 个 比特 的 水 印 
HE. 


e) 对 提取 出 的 一 维 水 印 数组 进行 升 维 、Arnold 逆 变 换 、 以 
图 像 格 式 输 出 ， 得 到 提取 出 的 水 印 图 像 。 


嵌入 水 印 的 JPEG 图 像 


获得 比特 解码 时 期 的 系数 Co~Ce 


Y 


Ace, SHE C +0 且 |C.|>1 的 系数 CG 


从 系数 的 最 低 有 效 位 中 提取 水 印 


Í N: 水 印 提取 完毕 


根据 统计 次 数 确定 当前 Y 
块 中 提取 的 水 印信 息 


| | 


选择 下 一 个 8*8 块 Arnold 反 变换 


Y 


获得 提取 出 的 水 印 图 像 


7 水 印 提取 流程 图 
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4 ”实验 结果 及 分 析 


41 实验 环境 


i5 处 里 器 ， 


本 文 测试 环境 为 : Windows10 Ji 
8GB 内 存 。 所 有 测试 均 在 


行 ,测试 过 程 中 为 了 便于 与 参考 文献 


使 用 了 


水 印 算法 中 经 典 的 载体 


图 像 ， 
Airplane、Boat， 大 小 均 为 512 X512， 


舰 版 操作 系统 ， 
Visual Studio 2012 下 运 
使 用 的 载体 图 像 进行 对 比 ， 
分 别 为 Lena, Pepper. 
量化 因子 为 100。 水 印 图 


Intel Core 


像 大 小 为 64X64， 原 始 水 印 图 像 如 


图 8 所 示 。 


经 过 多 次 实验 测 


WA, 大 的 值 取 5 时 效果 达到 最 佳 ， 


交 大 


图 8 原始 水 印 


实验 中 大 取 值 为 5。 


图 


Erud 


F1 
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不 可 感知 性 ， 
对 原始 图 像 载体 可 
方面 来 进行 评价 ， 
TURIN EE 


又 被 称 作 透明 性 ， 

见 性 的 影响 程度 ， 不 可 感知 性 一 般 通 过 两 个 
分 别 是 主观 和 客观 。 主 观 性 评价 是 
一 致 的 情况 下 ， 拘 入 水 印 后 的 图 


使 用 来 衡量 图 像 水 印 算法 


ul 


Xiv (F? 
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此 计算 ， 本 文 算 法 的 水 印 容量 为 0.078125bpp (bit per pixel) 。 
最 后 , 为 了 更 好 的 展现 本 算法 的 优势 , 本 文 与 文献 [10, 11] 进 行 
了 对 比 ， 对 比 结果 如 表 1 所 示 。 


表 1 本 文 算法 与 文献 算法 嵌入 容量 和 透明 性 对 比 
| 本 文 算法 文献 [10] 算 法 文献 [11] 算 法 
RARE/ RARE/ RARE/ 
512x512 PSNR PSNR 
10° bpp 10° bpp 10° bpp 

Lena 78.13 59.65 31.25 45.25 62.5 43.45 
Pepper 78.13 59.95 31.25 46.02 62.5 44.25 
Airplane 78.13 59.57 31.25 45.76 62.5 43.52 
Boat 78.13 56.34 31.25 44.93 62.5 45.06 


XA 1 可 以 看 出 ， 本 文 算法 在 答 入 容量 和 PSNR 上 均 远 高 
于 文献 [10，11] 的 算法 。 在 文献 [10] 中 ， 选 择 在 每 个 量化 后 的 块 
PRA 2 个 比特 的 水 印信 息 ， 由 于 直接 修改 量化 系数 后 在 哈 夫 


^E SS FS 
像 载 体 对 人 眼 


的 视觉 改变 量 越 小 代表 算法 的 透明 性 越 好 。 客 观 性 评价 一 般 通 
过 计算 峰值 信 噪 比 (PSNR ) 来 衡量 ,在 载体 图 像 和 由 入 容量 都 


一 致 的 情况 下 ， 原 始 载体 图 像 和 
PSNR 最 高 代表 水 印 算 法 的 透明 性 


嵌入 水 印 后 的 载体 图 像 的 
越 好 。 首 先 ， 从 主观 性 


上 进 


行 评价 , 如 图 9 和 10 所 示 , t 


体 图 像 和 


区 入 水 印 后 的 载体 


RZ, EARI SEEM 


无 任何 的 
图 像 的 痕迹 ， 说 明 算 法 的 透明 改 
上 看 ， 本 算法 对 于 Lena 图 像 ， 其 PSNR 为 59.6531; 
像 的 PSNR 为 59.9543; Airplane 图 像 的 PSNR 为 59.5651; Boat 
图 像 的 PSNR 为 56.3443， 可 以 看 
在 55 dB 以上, 而 同类 的 其 他 算法 , PSNR 基本 在 45 dB 左右 ， 


区 别 ， 没 有 任何 嵌入 过 水 印 


可 见 在 不 可 感知 性 方 


其 次 , 在 水 印 的 嵌入 容量 


E 良 好。 


其 次 ， 从 PSNR 的 结果 


Pepper 图 


H 4 张 测试 图 像 的 PSNR 都 


面 本文 算法 有 很 大 的 优势 。 


图 10 嵌入 水 印 后 的 载 


rh 


于 本 文 是 在 每 个 8X8 块 中 


选择 5 个 系数 来 嵌入 水 印 ， 所 以 每 个 块 中 会 嵌入 5 个 比特 。 


曼 编码 时 造成 了 误差 ， 因 此 PSNR 远 低 于 本 文 算法 ， 文献 [11] 
算法 利用 每 个 块 的 相关 性 来 生成 水 印信 息 ， 在 每 个 块 中 藤 入 4 


个 比特 ， 虽 然 文 献 [1H 在 藤 入 容量 上 与 本 文 差 距 不 大 ， 但 是 
PSNR 依然 远 低 于 本 文 算 法 。 另 外 ， 相 比 对 比 算法 来 说 ， 本 文 
算法 在 水 印 租 入 和 水 印 提取 更 加 简单 ， 同 时 由 于 不 需要 经 过 反 
量化 和 反 DCT 运算 即 可 提取 水 印 , 非常 适合 实时 应 用 。 提取 水 
印 运行 时 间 对 比如 表 2 所 示 : 
表 2 提取 水 印 运行 时 间 对 比 
本 文 算法 文献 [10] 算 法 文献 [11] 算 法 
载体 图 像 大 小 PS Bn m 
运行 时 间 /s 运行 时 间 /s 运行 时 间 /s 
256x256 0.608 1.806 1.986 
512x512 1.157 5.739 5.556 
1024x1024 1.155 13.716 16.894 
2048x2048 1.209 36.377 42.045 


UD 


其 中 ， 大 小 为 256X256 REE P BUB AE 32 X32 大 小 
水 印 图 像 , EAR UPS ELS RET UN. 64X64 大 小 水 印 图 像 。 本 文 
AEK EUN SUA AS EER, PERKEN EH RETI BUT 
解码 ， 并 且 一 旦 水 印 提取 完毕 即 可 停止 解码 过 程 ， 时 间 复 杂 度 
与 水 印 大 小 成 正 相 关 关 系 。 而 对 比 文献 算法 需要 进行 完全 解码 
后 才能 提取 水 印 ， 解 码 流程 复杂 ， 时 间 复 杂 度 高 
4.3 脆弱 性 测试 
对 于 脆弱 水 印 要 求 其 对 于 一 些 常见 的 处 理 
通过 提取 水 印信 


操作 比较 敏感 ， 
息 与 重 构 水 印 的 BER 和 NC 值 对 算法 的 脆弱 
为 测试 本 文 所 提出 方案 的 脆弱 性 ， 对 幅 入 水 印 后 
的 载体 图 像 分 别 做 了 如 下 几 种 攻击 操作 : 

a) 添加 方差 为 0.002 的 高 


性 进行 评价 。 


斯 噪声 ; 
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b) 添加 方差 为 0.002 的 椒盐 噪声 ， 也 肯定 不 同 ， 可 以 据 此 进行 窜改 定位 。 

c) 使 用 3X3 模板 的 低 通 滤波 ， 使 用 数组 wO 表示 原始 水 印 图 像 生成 的 一 维 水 印 数 组 ， 

d) 使 用 3X3 模板 的 中 值 滤波 ， mO 表示 从 嵌入 水 印 的 裁 体 图 像 中 提取 出 的 一 维 水 印 数组 ， 

O 重 压 缩 ， 量 化 因子 为 90; m) 为 原始 水 印 图 像 生成 的 一 维 水 印 数组 和 提取 出 的 一 维 水 印 

D 重 压缩 ， 量 化 因子 为 80 数组 的 差 值 并 取 绝 对 值 的 结果 ， 如 下 所 示 : 

图 11 给 出 了 本 文 算法 受到 以 上 攻击 之 后 的 提取 出 的 水 印 PEOR o) 
图 像 。 从 图 中 可 以 看 出 本 文 算法 受到 攻击 后 提取 出 的 水 印 图 像 。 其 中 ，i 代 表 图 像 中 的 第 i 块 ，m() 的 值 为 0， 那 么 说 明 丙 个 所 
均 被 破坏 ， 对 常规 攻击 都 比较 敏感 。 的 变化 不 大 ， 相 反 则 认为 该 图 像 英 已 经 发 生 改 变 ， 通 过 m) 即 

pom ss RE 可 以 定位 到 不 同 块 的 位 置 。 

d 为 了 测试 本 文 算法 针对 JPEG 图 像 被 恶意 审改 时 的 检测 和 


IRE ar ; ^d: 定位 方面 的 性 能 ， 本 文 进行 了 如 下 的 实验 操作 。 本 文 使 用 的 帘 
(a@) 攻 击 (1) 后 的 水 印 图 像 tcn en 图 像 〈c) 攻 击 (3) 后 水 印 医 改 图 像 均 为 Lena， 分 别 进行 了 3 次 窜改 测试 ， 第 一 次 测试 如 图 
3 L2) BER. E Lena 图 左下 角 用 Lena 的 头发 窜改 善 代 ， 第 二 次 
测试 如 图 12(b) 所 示 ， 帘 改 Lena 的 帽子 位 置 ， 第 三 次 测试 如 图 
12(c) 所 示 ， 分 别 窜改 Lena 的 左上 角 和 右 下 角 的 两 个 位 置 。 经 
过 帘 改 的 载体 图 像 相应 的 位 置 定位 如 图 13 所 示 ， 测 试 表明 本 


T EE PA di c in ro gr. 
(d) 攻 击 (4) 后 的 水 印 图 像 (e) 攻 击 (5) 后 水 印 图 像 ” (focii (6) KE É 


< 


3 
E E 文 算法 能 够 实现 对 审改 攻击 的 图 像 的 定位 。 


dé 3 给 出 了 经 过 上 述 攻 击 后 的 载体 图 像 提 取出 的 水 印 图 
像 与 原始 水 印 图 像 的 NC 和 BER 值 。 从 表 3 可 以 看 出 在 大 部 分 
攻击 下 ， 本 文 算法 NC 比较 小 ，BER 比较 大 ， 也 就 是 说 提取 出 
的 水 印 和 原始 的 水 印 差异 较 大 ， 可 见 本 文 算法 的 脆弱 性 较 好 。 
此 外 表 3 还 给 出 了 本 文 算法 与 文献 [10] 和 文献 [11] 算 法 的 对 比 
结果 ， 通 过 对 比 发 现 本 文 水 印 方案 在 受到 所 列 攻击 时 ， 绝 大 部 


JRE FARRA (b) 图 像 帽子 左上 被 窜改 (c) 图 像 左上 右 下 被 窜改 


(a 


图 12 Siu 


分 攻击 操作 下 的 NC 均 小 于 文献 [10,11]，BER 均 大 于 文献 
[10,11]， 说 明 本 文 所 提 的 方案 的 脆弱 性 要 优 于 对 比 文献 。 
表 3 本 文 算法 和 文献 算法 脆弱 性 对 比 
攻击 类 型 本 文 算法 文献 [10] 算 法 文献 [11] 算 法 


BER NC BER NC BER NC 
图 13 帘 改 定位 检测 结 


无 攻击 0 1 0 1 0 1 

高 斯 噪声 0.5090 0.5853 0.2021 0.7905 0.4823 0.5923 5 ”结束 语 

椒盐 噪声 0.1609 0.8811 0.4406 0.3236 0.4175 0.4750 

低 通 滤波 0.5012 0.5886 0.4266 0.5818 0.4703 0.6270 AK XCRU EDS LER BOR ES EUIS IER I T CERA UNE 


性 ， 提 出 了 一 种 基于 炳 编码 的 JPEG 压缩 域 脆弱 水 印 算法 。 本 
算法 选择 将 水 印信 息 元 余 嵌 入 到 每 个 图 像 块 的 个 DCT 交流 
系数 炉 编 码 后 的 最 低 有 效 位 中 ,一 方面 可 以 提高 算法 的 稳定 性 ， 
另 一 方面 还 可 避免 量化 损失 和 哈 夫 曼 编码 误差 等 问题 ， 使 得 算 
法 透明 性 较 高 。 实 验 结果 表明 ， 算 法 具有 较 大 的 嵌入 容量 和 很 

好 的 不 可 感知 性 ， 对 常规 攻击 脆弱 性 较 高 且 能 够 对 恶意 窜改 进 
时 期 一 元 组 中 右边 的 参数 来 进行 水 印 提取 ， 如 果 载 体 图 像 的 内 行 准确 定位 。 算 法 实现 过 程 简单 ， 仅 需 对 编 解 器 进行 较 少 改动 
9 窜改 ， 那 么 经 过 解码 后 获得 的 相应 参数 值 不 即 可 实现 ， 在 实现 水 印 的 盲 提 取 的 同时 保证 了 算法 的 效率 ， 提 
会 改变 。 如 果 载 体 图 像 被 恶意 窜改 ， 那 么 窜改 区 域 图 像 块 中 相 。 ”高 了 算法 的 实用 性 
应 获得 的 参数 会 大 幅 的 改变 ， 从 而 从 该 窜改 区 域 提取 出 的 水 印 


中 值 滤波 0.5154 0.5702 0.5083 0.5969 0.4056 0.5855 


重量 化 90 — 0.5471 0.5172 0.1064 0.8896 0.1186 — 0.8513 
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0.5425 0.4987  Á 0.2021 0.7886 0.1895 0.9432 
4.4 窜改 定位 测试 
本 文 算法 是 基于 哈 夫 曼 解 码 和 行程 解码 后 , 获取 比特 编码 
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